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We investigate the real-time dynamics of the chiral magnetic effect in quantum electrodynamics (QED) 
and quantum chromodynamics (QCD). We consider a field configuration of parallel (chromo)electric and 
(chromo)magnetic fields with a weak perpendicular electromagnetic magnetic field. The chiral magnetic effect 
induces an electromagnetic current along this perpendicular magnetic field, which we will compute using linear 
response theory. We discuss specific results for a homogeneous sudden switch-on and a pulsed (chromo)electric 
field in a static and homogeneous (chromo)magnetic field. Our methodology can be easily extended to more 
general situations. The results are useful for investigating the chiral magnetic effect with heavy ion collisions and 
with lasers that create strong electromagnetic fields. As a side result we obtain the rate of chirality production 
for massive fermions in parallel electric and magnetic fields that are static and homogeneous. 



I. INTRODUCTION 

In quantum electrodynamics (QED) and quantum chromo- 
dynamics (QCD) certain gauge field configurations will in- 
duce chirality as follows from the axial Ward identity 0][2). 
These gauge field configurations have a net topological charge 
density, which means that the contraction of the field strength 
tensor with its dual is nonvanishing. An example of such 
gauge field configuration is one in which the (chromo)electric 
and (chromo)magnetic fields are parallel. For massless parti- 
cles, chirality is the asymmetry between the number of parti- 
cles plus antiparticles with right-handed helicity and the num- 
ber of particles plus antiparticles with left-handed helicity. 
Here right-handed helicity means that spin and momentum 
are parallel, whereas left-handed helicity means they are anti- 
parallel. If a system of charged fermions possesses a net chi- 
rality, then an applied (electromagnetic) magnetic field will 
induce an electromagnetic current in the direction of the mag- 
netic field [3j. This is called the chiral magnetic effect. In 
magnetic fields that are large enough to fully polarize the 
fermions, one can easily convince oneself that the magni- 
tude of this current is equal to the chirality times the absolute 
charge of the fermion J4] . 

To improve the understanding of the chiral magnetic effect 
it is important to obtain the magnitude of the induced elec- 
tromagnetic current in different situations. In many studies of 
the chiral magnetic effect the chirality is introduced by hand 
through a chiral chemical potential /15 [3] [5] |6] . However, in 
order to get a full understanding of the chiral magnetic effect, 
one should include the dynamics that leads to a net chiral- 
ity. Such dynamical studies have been performed in Euclidean 
space-time using lattice QCD [7| and with analytic methods, 
both in Euclidean [8 9| and Minkowski [ 1 ] space-time. 

In this article we will study the real-time dynamics of the 
chiral magnetic effect with an analytic method. We will con- 
sider a specific gauge field configuration of parallel homoge- 
neous (chromo)electric and (chromo)magnetic fields that are 
pointing in the z-direction. These parallel fields are the source 
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FIG. 1: Field configuration in which the chiral magnetic effect is 
studied in this article. The electric field along the z-direction (E z ) is 
homogeneous and can have arbitrary time dependence. The magnetic 
field along the z-direction is static and homogeneous. The goal is 
to compute the current density j y along the magnetic field in the y- 
direction (6 y ). 



of the chirality. The (chromo)magnetic field is assumed to be 
constant in time, but the (chromo)electric field can have any 
time dependence. Perpendicular to these fields we will con- 
sider an electromagnetic magnetic field with arbitrary time de- 
pendence that points in the y-direction. We have sketched the 
situation in Fig. [T] The goal of this article is to obtain the 
current density in the y-direction, which is due to the chiral 
magnetic effect. 

The chiral magnetic effect will take place unsuppressed if 
the fields have a magnitude that is at least of order of the 
fermion mass squared. In the near future, one hopes to cre- 
ate such enormous electromagnetic fields with high intensity 
lasers (see e.g. ifTTI ). If one could engineer a field configura- 
tion like the one of Fig. [T]with these lasers, one could falsify 
the chiral magnetic effect in QED in a controlled environment. 
In heavy ion collisions a similar configuration built out of flux- 
tubes of parallel chromoelectric and chromomagnetic fields is 
induced [12]. Perpendicular to these color fields a large elec- 
tromagnetic magnetic field is generated through the currents 
created by the colliding charged ions (H [13] . 

In QCD the probability of generating gauge field configu- 
rations with positive and negative topological charge is equal 
(assuming = 0). Hence in heavy ion collisions, the cur- 
rent along the perpendicular field vanishes on average. How- 
ever, because of fluctuations in the initial state lfl2l . and the 
presence of sphalerons [14| the topological charge fluctuates. 
Hence the current along the perpendicular magnetic field will 
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fluctuate as well. That will lead to differences in net electric 
charge between the two sides of the reaction plane (the plane 
in which the beam-axis and the impact parameter lie) on an 
event-by-event basis 1 15 1. This effect can be analyzed experi- 
mentally through a charge correlation study lfl6l . Charge cor- 
relations have indeed been observed at both the Relativistic 
Heavy Ion Collider (RHIC) 02) and the Large Hadron Col- 
lider (LHC) IfTSl . However, currently it is not clear if these 
correlations are caused by the chiral magnetic effect or by an- 
other mechanism lfT9l . A possible observation of the chiral 
magnetic effect in QCD can teach us about the relevance of 
configurations with topological charge. 

Theoretically, the field configuration displayed in Fig. [T] is 
interesting, because when all fields are static and homoge- 
neous it is possible to obtain a full analytic expression for 
the induced current density along the perpendicular magnetic 
field IflOl . The current density behaves exactly according to 
the qualitative expectations of the chiral magnetic effect. For 
example, for a large perpendicular magnetic field the magni- 
tude of the current was found to be equal to the chirality times 
the absolute value of the charge of the fermion. Furthermore, 
the current density is suppressed by the fermion mass, and 
vanishes if there is no chirality production. 

The calculation in Ref. ifTUl was build upon a nice argument 
based on Lorentz transformations and the particle production 
rate in parallel, homogeneous and static fields. However, this 
argument cannot be applied to more general situations that are 
useful to investigate. For example, in heavy ion collisions the 
magnetic field has a strong time dependence. Furthermore, 
with high-intensity lasers one will create inhomogeneous and 
time-dependent fields. And for purely academic reasons it is 
of interest to analyze the real-time dynamics of the chiral mag- 
netic effect in a spherical symmetric background gluon field 
that carries one unit of topological charge, and is a solution of 
the Yang-Mills equations. 

To obtain the induced current density along the magnetic 
field in a general situation one can in principle just solve the 
Dirac equation numerically, and construct the current from the 
solutions. However, this might be quite difficult in practice. In 
order to make progress, we will start a program to investigate 
the chiral magnetic effect using linear response theory. We 
will expand to first order in the electromagnetic field that is the 
source of the perpendicular magnetic field. In this way we can 
investigate the chiral magnetic effect for weak perpendicular 
magnetic fields that have arbitrary space and time dependence. 
In practical terms we will confirm in this article the results of 
Ref. IflOll in the weak magnetic field limit. Furthermore we 
will analyze the chiral magnetic effect in pulsed electric fields. 

Most of the calculations in this article will be performed 
in QED. From the QED results it is easy to obtain the results 
for QCD as in Ref. [ 10 1 and as explained at the end of the 
introduction. We hope that our results will inspire the reader 
to perform other studies of the chiral magnetic effect using 
linear response theory. 

This article is organized as follows. In Sec.[II]we will con- 
struct the fermion propagator of QED in a background of par- 
allel electric and magnetic fields. To support Ansatze made 
in Ref. IfTOl . we compute the electromagnetic current density 



for massive fermions analytically in Sec. Ill In this section 



we will also obtain an analytic expression for the rate of chi- 
rality production for massive fermions in parallel electric and 



magnetic fields. In Sec. IV we will then apply linear response 



theory to obtain the induced current density along a magnetic 
field that is perpendicular to these background fields. We will 
discuss our numerical results in Sec.|V]and show among other 
things that we can reproduce the results of Ref. [10] in the 
limit of weak magnetic fields. The conclusions of this work 
can be found in Sec.lVIl 

Some details of this work can be found in the appendices. In 
Appendix[A]we give several relations involving the wavefunc- 
tions in a magnetic field. Finally in Appendix [B] we discuss 
the wavefunctions in a sudden switch-on electric field. 

The notational conventions we will use in this article are 
summarized below. We also explain the relation between the 
QED and the QCD setup below. 



Notation 

We will use as a metric g^ = diag(l, -1, -1, -1). We write 
p — p>* — (p°, p) to denote a 4- vector. A component of a four- 
vector will be written with either numbers or with an italic 
sub- or superscript. For example, p z = p 3 = -p z = -p$. 
We will use roman subscripts to indicate a component of a 3- 
vector, i.e. p = (p x ,p y ,p z ) = (p x ,p y ,p z ). Hence p z = p z = 



-p z . Likewise we write x 
simplify to x = (x, y, z). 



(x x , x y , x z ) which we most often 



Relation between QCD and QED setup 

As in Ref. 1 10 1 we consider an Abelianized gluon field con- 
figuration of the form A° = in the QCD setup. Here 
n" denotes a direction in the adjoint space. Since the elec- 
tromagnetic current density, the object we wish to compute, 
is gauge invariant, we can perform a gauge transformation to 
rotate if to the a — 3 direction. From the covariant deriva- 
tive Dfj - df! - ig&nt 3 , it then follows that the red quarks 
effectively experience an electromagnetic field of magnitude 
qA^ - —gtRplI whereas the green quarks experience an effec- 
tive electromagnetic field with size qA^ - gJ^^/2. The blue 
quarks are not affected by the field. Here q is the charge of the 
fermion and g is the QCD coupling constant. The current den- 
sity along the perpendicular magnetic field in the QCD setup 
can thus be obtained from the QED result, yielding 

jy.QCD = 2 2_j jy,QED(qE z = g& z /2, qB z = gS z /2, 

f 

q = qf,m = nif). (1) 

Here / is a sum over quark flavors, m is the mass, and 
{& Z )E Z and {T> Z )B Z denote respectively the z-component of the 
(chromo)electric and (chromo)magnetic field. 

There is an important difference between the QED and 
QCD setup. In QED the electromagnetic current has also a 
component in the z-direction due to acceleration of charges 
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along the electric field. In the QCD setup however, this contri- 
bution to the electromagnetic current vanishes due to a cancel- 
lation from the contribution of red and green quarks. The elec- 
tromagnetic current in the QCD setup has thus only a compo- 
nent in the y-direction due to the chiral magnetic effect. 



II. FERMION PROPAGATOR IN ELECTROMAGNETIC 
BACKGROUND 

The goal of this section is to obtain a convenient expres- 
sion for the fermion propagator in a background of electric 
and magnetic fields both pointing in the z-direction. We will 
consider a homogeneous static magnetic field, B = B z e z . The 
electric field is assumed to be homogeneous as well with an 
arbitrary time-dependence, E = E(t)e z . The corresponding 
gauge fields will be chosen as: A — and A = (0, B z x, A z (t)) 
with lim,^_ M A z (t) = so that E(t) = -d,A z (t). 

In this section we will first compute the Dirac spinors in 
this electromagnetic background. Using these spinors we will 
construct the quantum field for fermions, from which we will 
derive the different fermion propagators. 



A. Dirac spinors 

The spinors tf/(x) can be obtained from the Dirac equation 
which reads 



(iyDf, - m) i//(x) = 0, 



(2) 



where the covariant derivative - 8^+ iqA^ with q the elec- 
tric charge of the fermion. 

To solve the Dirac equation in the background given above 
we will apply the Ritus method [20]. This method amounts to 
making the following ansatz for the particle spinors (+) with 
momentum p and antiparticle spinors (-) with momentum —p 

tf, ps (x) = F + p (t)G p (xyw +1 ?< z u s (p + ), (3) 
if, ps (x) = F-(t)G p (x)j™ +ip *v s {p-). (4) 

Here the four-vector p 4 ^ = (k, 0, +A, 0), and u s (p+) and v s (/L) 
denote the usual particle and antiparticle spinors of the free 
Dirac equation with momentum p ± . The free Dirac spinors 
satisfy the following equations 



m)u s (p+) = 0, 



+ m)v s (pS) — 0. 



(5) 



Furthermore, k and A are constants that depend on p, and will 
be determined below. Finally, F p (t) is a 4 x 4 matrix that has 
to commute with y 1 and y 2 , and G p (x) a 4 x 4 matrix that 
should commute with y° and y 3 . 

By inserting the Ritus ansatz into the Dirac equation it fol- 
lows that this ansatz is indeed a general solution if the matrices 
F p (t) and G p (x) obey the following two equations 

[iy () d t - 7 3 ( Pz -qA z (t))]F;(t) = ±KF p (t)y°, (6) 



We will solve these two equations explicitly below. 

By multiplying Eq. |5]l from the left with (j> ± + m) it fol- 
lows that p\ = m 2 , so that k 2 = A 2 + m 2 . We will normalize 
the free Dirac spinors as in Ref. ED . In that convention the 
polarization sums, which are the only non-trivial property of 
the free spinors required in this article, read 



2^jU s {p)u s {p) = {f + m), 

s=± 



(8) 
(9) 



In the remaining part of this subsection, we will determine 
the matrices F p ( f) and G p (x), and the values of k and A ex- 
plicitly. Without loss of generality both F p (t) and G p (x) can 
be taken diagonal. 



1. Computation qfG p (x) 

From Eq. (|7ji it follows that G p (x) is a function of x - 
p y /qB z . Therefore let us write the diagonal elements of G p (x) 
as g a (x - p y /qB z ) with a - 1 . . .4. By multiplying Eq. (|7]) 
from the left with y 2 one finds that 



(3 X + qB z x)gi(x) 
(<9 X -qB z x)g 4 (x) 



^4(x), 
-Agi(x). 



(10) 
(11) 



Since the other two components satisfy the same set of equa- 
tions, g2(x) is proportional to g4(x) and gi(x) is proportional to 
gi(x). In order to ensure that G p (x) commutes with y° and y 3 
as required, we have to take giix) = g4(x) and gi(x) = g3(x). 



By combining Eqs. ( 10 1 and (Hi one finds that g a (x) satis 



fies the following equation 

(-d 2 x + q 2 B 2 z x 2 - sqB z )g a (x) = A 2 g a (x), (12) 

with s = 1 for a = 1 , 3 and s = - 1 for a = 2, 4. Equation 
( p"2| ) is the eigenvalue equation of a harmonic oscillator with 
angular frequency \qB z \. The eigenvalues are given by A 2 

(In + l)\qB z \ - sqB z with quantum number n = 0, 1, 2, The 

corresponding eigenfunctions (f>„(x) read 



(p„(x) 



1 



V2"n! 



\qB z 



1/4 



H n {^I\q^\x)tx^L l -\qB z \x 2 \ 



(13) 

here H„(z) denotes a Hermite polynomial of degree n. For 
convenience we define <^_i(x) = 0. The eigenfunctions <p n (x) 
are normalized as 



x: 



dx(p„(x)(p m (x) = 8„ m , 
and satisfy the completeness relation 

CO 

2^ 0„(x)0„(x') = S(x - x'). 



(14) 



(15) 



re=0 



Let us now write A 2 = 2k\qB z \ with k = 0, 1,2, The full 



[iy l d x -y 2 (p y -qB z x)]G p (x) = AG p (x)y 2 . (7) solution to Eq. (|12[i can be found by using that g\ (x) and g^(x) 
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must have the same value of A 2 . This gives for qB z > 0, 
gi(x) = ci$fc(x) and g4(x) = C40jt_i(x). Choosing the nor- 
malization such that ci = 1, it follows from Eq. ( fT0) > that 
C4 = c\ = 1. In this way one obtains also the sign of /I, giving 
A = sgn(gB z ) ^j2\qB z \k. As long as gB z > one then finds 
ga(*) = <At(x) for a = 1,3 and g fl (x) = <At-i(x) for a = 2,4. 
If qB z < the solution is g n (x) = ^-[(x) for a = 1,3 and 
g a (x) = ^(x)fora = 2,4. 

The diagonal matrix G(x) can now be written in the follow- 
ing compact notation 



(16) 



where we introduced g p+ (x) 
0*-i( x 

^ = (l+i5Sgn(^B z W)/2. 



^ x ~~ w) and ^p-W = 
■^r). The projection operators P s g are given by 



2. Computation of F^(t) 

Let us write the diagonal elements of the matrix Fp(t) as 
/^(f) with a = 1 ... 4. Multiplying Eq. (|6]) from the left with 
y° yields 



[id t +p z -qA z (t)]Jf(t) 
[id t - Pz + qA z (t)]Jf(t) 



(17) 
(18) 



Since the other two components satisfy the same set of equa- 
tions it follows that f*(t) oc and oc /*(?). In order 
to ensure that commutes with y 1 and y 2 we have to take 
/±(f) = Jf(t) and /*(?) = /f (t). We can therefore write 



(19) 



where /*+(?) = jf (i), fpAO = f^(t) and the projection opera- 
tors 9) = (i + sr 3 r°)/2. 



By taking the complex conjugate of Eqs. ( 17 1 and ( 18 1 it 
can be seen that f ps (t) satisfies the same two differential equa- 
tions as /p_ s (0*- Hence both functions are proportional to 
each other. Since it is natural to normalize the particle and 
antiparticle spinors in the same way, we take 



(20) 



We will use this relation throughout to express our final results 
in terms of particle wavefunctions only. 
From Eqs. ( fTTj ) 
_i_ i f± tt\d 



and 



18 1 it also follows that the sum 
is independent of time. As it turns out it 



l/| + (0l 2 + l/±_(0l 

will be convenient to normalize this combination as 

\fp + (t)\ 2 + l/p-(0| 2 = 2. 



(21) 



In the next subsection we will check that the normalizations 
we have made are consistent, by verifying that the quantum 
fields satisfy the canonical anti-commutation relations. 

In general the wavefunctions f ps (t) can only be obtained 
numerically. There are only a few cases in which analytic 
solutions are known. We will now discuss a few of them. 



For vanishing electromagnetic field, A z (t) - 0, combining 
Eqs. ( fT7| ) and ( fT8| ) gives 



[-%-pl]f P t s (t) = '?f P t s (t), 



(22) 



with s = +. The solution of the last equation is a linear combi- 
nation of phase factors exp(±ipoO, where po = -\/k 2 + p 2 . The 
two different solutions correspond to particles and antiparti- 
cles, hence ff (f) = exp(+ipoO where is a normalization 



constant. From Eq. ( 17 1 it follows that the ratio of the normal- 
ization constants is c%/ct = k/(pq ± p z ). Applying Eq. (21 
gives (c*) 2 + (c*) 2 = 2, so that 



Po + sp z 



Po 



exp(+i/? f)- 



(23) 



If k = 0, the two differential equations for fi s (f) decouple. 



In that case Eqs. 
wardly, yielding 



17 1 and ( 18 i can be integrated straightfor- 



fp t s®\ K =o = ^ ex P ^W-i* f dt'qA z (t') 

*J —CO 



6(+sp z ). 



(24) 

In a sudden switch-on electric field the wavefunctions f ps {t) 
are known analytically J22) . We review the calculation in Ap- 
pendix[B] It is also possible to obtain an analytic solution in a 
pulsed field of the form E z (t) = E z / cosh 2 (f/r) 11231 . 



B. Quantum field 



The Dirac field in the electromagnetic background is given 

by 

s=± p V P 

where we introduced £ p = Yi7=o f f IF anc l K p ~ K - Here 
i/'p^jc) are the particle (+) and antiparticle (-) spinors in the 
background field, which are given explicitly in Eqs. ([3]) and 
Q. The operators b ps and d ps denote respectively the annihi- 
lation operators for particles and antiparticles with momentum 
p and spin s in a background magnetic field. The creation op- 
erator for antiparticles in Eq. ( |25] l has negative momentum, 
reflecting the fact that in our notation (A ps (x) denotes an an- 
tiparticle spinor with momentum -p. 

The creation and annihilation operators satisfy the follow- 
ing anti-commutation relations 

\b , bt ,} = {d , ct , ,} 

1 ps ' p's" ' ps ' p's" 

= (2n) 2 5 kk >5{p y - p' y )6(p z - p' z )8 ss >. (26) 

All other anti-commutation relations vanish. 

To check that all normalization conditions are consistent we 



will verify that the quantum field given in Eq. ( 25 i satisfies the 
canonical equal-time anti-commutation relation, which reads 



{W a (t,x)M(t,x')} = 5 ab 6 3 (x-x'). 



(27) 
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Using the explicit expression of the quantum field gives after 
using the properties of the creation and annihilation operators 

mt, X ),w\t,x')}= 2Zr~^ (f,a! W'' a! ' )t - (28) 



u,s=± p 



Inserting the explicit solution for the spinors, and summing 
over spins by applying Eqs. ([8]l and (|9|l yields 



w=± p 



x [ U F u p (t) 7 °F u p (tf [mG p (x)G p (x') - AG p (x)y 2 G p (x')] 

+ K p F l p (t)F p (t^G p (x)G p (x')}. (29) 

The last equation can be simplified by inserting the explicit 
expression for F p (t). In this way we find that 

\j]F p (i)F p (rf = U, (30) 

u=± 

YjUF^trfFpit? =0. (31) 

u=± 

To obtain the last two equations, we have used that by com- 



bining Eqs. (20 1 and (21 1 it can be shown that 



l/ P + s (0| 2 + \f p M 2 = 2, 



(32) 



&v)fe,(ty = f ps (t)f P s(tr ■ (33) 

Inserting Eqs. ( |30] > and pT[ ) into Eq. (29 1 yields 

{¥(*, x), ¥'(t, x')} = Yj e i ' ,y(y_y ' )+ift(z_z ' ) Gp(x)G p (x'). (34) 



By using the explicit expression for G p (x) and applying the 
completeness relation Eq. (15 1, the canonical anticommuta- 



tion relation, Eq. ( 27 1, follows directly. Hence the normaliza- 



tions we have chosen are consistent. 



C. Propagator 



By inserting the explicit expressions for the Dirac spinors, and 
summing over spins, the two-point correlation functions be- 
come 



S ± (x, x 



a i/>y(y-y')+ip 2 (z-z') 



x F;(t)G p (x)(p ± ± m)y F±(t'y 7 °G p (x'). (38) 



To evaluate the current density and related quantities, one 
has to contract a two-point correlation function with an arbi- 
trary combination of gamma matrices denoted by F. These 
quantities can be expressed in terms of S ± (x, x') in the fol- 
lowing charge symmetric way, 

<0|¥(f,sc')ri'(f,x)|0> = ~Y j utx[S u (t,x;t,x')r] 

u=± 

+ ^tx[y°T]6(x-x'), (39) 

For the current density f we have Y = qyf, for the chirality 
«5 we have F = y°y 5 , and for the pseudoscalar condensate 
F = iy 5 . Only for F = y° the second trace in Eq. (39 1 does not 
vanish. 



III. CURRENT AND AXIAL ANOMALY IN PARALLEL 
ELECTRIC AND MAGNETIC FIELD 

We will now compute the induced current density, chiral- 
ity density and pseudoscalar condensate in the background of 
parallel homogeneous time-dependent electric and static mag- 
netic fields. To compute these quantities we need to evaluate 
the trace of S ± (x, x')T as follows from Eq. ( 39 1. This trace can 



be easily performed using a symbolic manipulation program 
such as Mathematica. Inserting the explicit forms of F p (t) and 
G p (x) gives 

trfs^JCxV] = 4Z Z e iMy - y ' )+ipAz - z,) 

p r,s=± 

x sf; s {t)f; s {ty gpr (x) gpr {x'), (40) 



Let us introduce the following definitions for the two-point 
correlation functions 



S + ah (x,x') = (0\V a (x)%tf)\0), 
S- ah {x,x') = {Q\V b {x'yV a {x)\0). 



(35) 
(36) 



Here |0) denotes the in-vacuum, which in this article is the 
vacuum before the electric field has been switched on. The 
different propagators (retarded, advanced, Feynman) can be 
found by the appropriate linear combinations of the two-point 



functions. By applying Eq. (25 i it follows that the two-point 



correlation functions expressed in terms of Dirac spinors read 
S± ab ^ *') = Z Z aT W x '^y\ ■ (37) 

s=± p P 



trfS^JCxVy 5 ] = Z Z e^ (y - y ' )+ift(z - z,) 

p r, s=± 

X s^{qB,)sf^{t)f; s {t'y rgpr (x)g pr {x'), (41) 



tr [£*(*, jc')7 5 1 = +Z Z ^e^ (y - y '> +i ^ z - z '> 

p r,s=± 1 

x sgn(qB z )sf^t)f p -Jtrrg pr (x)g pr (x). (42) 

We have to evaluate these correlators at equal time (f = t') 
with the x and y-components of the direction vectors equal 
(x = x' and y = y'). For reasons we explain below, we will 
keep z and z' different and introduce A = z- z'. 
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In this situation we can perform the p y integration using the 
orthogonality properties of the functions g ps (x). This yields 

-^«tr[r(f,A) r 3 ] = 



An 



\J n _|_ 



-i^«tr[5"(f,A)r 5 ] = 

u=± 

sgn(m) S f£ el " A Z<«/p + -^*L- (45) 

where a* = 1 for k = and = 2 for > 0. We have 



made use of Eqs. ( |20| > and (21 1 to express all results in terms 
of particle wavefunctions f£ s (t). 

The quantities we will compute can be obtained from 
Eqs. (f43)-(|45]> in the limit A -> 0. Naively putting A = will 
not give a well defined result, due to the presence of ultraviolet 
divergences. Therefore we have to regularize our expressions. 
For consistency, this regularization should be performed in a 
gauge invariant way. 

A natural way to achieve this is by using the point-split reg- 
ularization |24|. Instead of evaluating the two-point functions 
at z = z' one computes them at z - z' = A and integrates them 
over a distribution /i(A). This distribution has to be normal- 
ized to unity, and should be sharply peaked around A = 0. We 
will choose h{A) = exp(-A 2 /4e)/(2 -\fne) and take the limit 
e -> 0. 

In order to maintain gauge invariance, the correlators 
5 ± (f, A) have to be augmented by a gauge link link U con- 
necting z with z' ■ This gauge link is for both S + (t, A) and 
S~(t, A) given by 



U(t, A) = exp 



iq I dx?A M (x) 



exp [-iqA z (t)A] . (46) 



Summarizing, the full point-split regularization prescription 
reads 

tr, [S ± (0T] = lim j dA/i(A)t/(f, A) tr [S ± (r; A)T] , (47) 

here the subscript r stands for regularized. 

We will now apply this regularization prescription to 



Eqs. (43 H 45 1 to evaluate the current density, chirality den- 
sity and pseudoscalar condensate. The integral over A can be 
performed exactly giving 



\qB z \ 
'i-^'6' 1 2n 



fe=0 



x [l/ p + + (f)| 2 - l] , (48) 



"5(0 



lim** 

mo In 



I 



d P z p-e[fe-gA.(t)] 2 
In 



X [l/; + «lLo - l] . (49) 



qB z 

77(f) - lim i sgn(m) 

e->0 4n 



I 



In ' 



■ 2 ^(0^(0*1^, (50) 



where we used that from Eq. (21 1 it follows that £. s |/p S (f)l 2 



2(\fp + (t)\ - 1)- In the absence of a regulator and in the limit 
of vanishing magnetic field, Eq. ( |4"8j ) agrees with the results 
obtained in Ref. |25|. 

We can now verify the axial anomaly relation. In the mass- 
less limit one can show that \f£+(t)\l =0 - 1 = -sgn(/? z ). We 



can now perform the p z integral in Eq. (49 1 in the limit e 
giving 







"5(0 = -^B z A z (t). 



(51) 



Taking the derivative with respect to time in the last equation 
gives the axial anomaly relation for massless particles in par- 
allel electric and magnetic fields 



dris(t) q 



2n 2 



B z E z (t). 



(52) 



If the fermions are massive the anomaly relation contains an 
additional term proportional to the pseudoscalar condensate. 
To obtain the anomaly relation in this case we will perform 
the time derivative on the chirality given in Eq. ( |49| ) explicitly. 
The limiting procedure e — > is equivalent to performing 
the integration over p z in an interval symmetric around p z = 
qA z (t). Therefore 



f A v vb z r^' AA,) d Pz 

n 5 (f) = lim | 

a^oo In J-A +qAl ( t) 2n 



[l^WlLo-l]. (53) 



The time derivative of n 5 (f) contains a part arising from the 
integration boundaries and a part from the derivative on the 



wave-functions. It follows directly from Eq. ( 17 1 that 



d,\f P+ (t)\ 2 = ik p 2 s/; s (t)f;_ s (ty. (54) 



Also we will use that k p \^=o = \m\. Furthermore, for large p z 



one can neglect qA z (t), and from Eq. ( 23 i it can be shown that 
lim^oc \f p +(t)\ 2 = and lim^-M \f p +(t)\ 2 - 2. By using 
Eq. (50 1 and applying the time derivative to the integration 



boundaries it follows that 



dn 5 (f) 
dt 



2mr]{i) + ^B z E z {i). 



(55) 



which is exactly the anomaly relation in the presence of mass. 

Let us now consider a sudden-switch on electric field of the 
form E(t) = E z 6(t). We discuss the functions f ps {t) in this 
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field in Appendix [B] We can only compute the induced cur- 
rent density and the chirality density analytically in the large t 



limit. Applying Eq. (B13 1 we find that for large t 



dm 

dt 



dn 5 (t) 
dt 



1 



\qB z \qE z __^ ^ _-2wfc|g| 



In 1 



I 



k=o 

l B, 



In 1 



(56) 
(57) 



Since these results are obtained for large t in a sudden-switch 
on electric field, they are exact in a constant electric field. 



The result Eq. ( 56 1 was derived analytically in a different 
way in Ref. [29 1. It was also proved to be correct numerically 
in Ref. Il30l . The result is easy to understand starting from 
the production rate of fermion antifermion pairs in parallel 
homogeneous electric and magnetic fields HI OH . That rate per 
unit volume equals [22, 31] (see also [32-34)) 



r = 



q 2 E z B, 
4n 2 



For B z — this rate reduces to the pair production rate in a ho- 
mogeneous electric field, that was first obtained by Schwinger 
HI. The production of pairs gives rise to a homogeneous cur- 
rent density that has to point in the z-direction because of sym- 
metry reasons. The particles are accelerated continuously by 
the electric field. Therefore, at some point they will reach (al- 
most) the speed of light. Hence, every time a pair is created 
the current will eventually grow by twice the charge of the 
fermion. So therefore the rate of change of the current density 



is given by d t j z = 2qYsgn(qE z ). Inserting Eq. (58 i we see that 



we exactly recover Eq. (56i. One can also use this argument 
in the opposite order, in order to derive the pair production 
rate from the calculation of the current density. 

The result Eq. ( |57] i generalizes the well-known production 
rate of chirality in parallel electric and magnetic fields for 
massless fermions to massive fermions. We are unaware of 
an earlier derivation of this result. The mass suppresses the 
production of chirality. By combining Eq. ( |5"7| ) with Eq. ( |55| l 
we find that the pseudoscalar condensate in static, homoge- 
nous and parallel electric and magnetic fields equals 



7(0 = 



q 2 E z B z 
An 2 m 



^e » - lj. 



(59) 



IV. LINEAR RESPONSE TO MAGNETIC FIELD 

In the previous section we have considered a time- 
dependent electric field and a constant magnetic field that 
were both pointing in the z-direction. To this field config- 
uration we will now add a time-dependent magnetic field in 
the y-direction, denoted by B y (t). This magnetic field will be 
accompanied by a perpendicular electric field as can be seen 
from Faraday's law, Vx£ = -dB(t)/dt. The additional mag- 
netic field will induce a current density in the y-direction. In 



this section we will compute this current density to first order 
in B y (t) using linear response theory. In the next section we 
will use this result to study the chiral magnetic effect. 

Let us write the full electromagnetic field as A^{x) = 
A M (x) + A /J (x). Here A^{x) denotes the background field, con- 
sisting of the electric and magnetic fields pointing in the z- 
direction. The field A^(x) denotes the perturbation on this 
background, which in this case is the magnetic field in the y- 
direction with its corresponding perpendicular electric field. 

From linear response theory it follows that to first order in 
A^ix) the induced current density in the electromagnetic field 
Af(x) equals f(x) = %(x) = f A (x) + Sfi^x) where 



= / d A x'n^{x,x')A v {x'). 



(60) 



Here the retarded current-current correlator (or equivalently 
photon polarization tensor) in the background field A^x) is 
given by Tl^(x, x') = U^(x, x')6(t - t') with 

Yl^{x,x) = -iq 2 (0\ TGOy^Oe')] |0>. (61) 



Using Eqs. (25 1 and 07b we can express TVt (x, x') as 



Yl^(x,x') = -\q 2 Y J uK[fS u (x,x)y v S- u (x ,x)} , (62) 

«=± 

where S ± (x,x') is the two-point correlation function in the 



background field, given explicitly in Eq. (38 1. 

Since the background electric and magnetic fields are both 
pointing in the z-direction, they cannot solely induce a current 
in the y-direction. As a result f^(t) = 0. Hence the induced 
current density in the y-direction can only arise from the per- 
turbation and is therefore of the following form 



MO 



•J — oc 



dt'H(t, t')By(t'). 



(63) 



where H(t, f) can be obtained from Tl^ix, x') as we will ex- 
plain in more detail below. The photon polarization tensor 
in an electric plus magnetic background has been studied by 
other authors before in different contexts (26). Furthermore, 
the photon polarization tensor in a purely magnetic back- 
ground has been studied in detail in several works 11271 . for 
recent analyses and applications we refer to I5l l28l . 

In the following subsections we will compute the func- 
tion H(t, f) in two cases, labeled by A and B. In case A we 
will take the only non-vanishing component of the perturba- 
tion field to be A z (x) = -B y (t)x. In case B the only non- 
vanishing component is chosen as A x (x) = B y (t)z. These 
two cases lead to the same magnetic field in the y-direction, 
B y (t), but give rise to different perpendicular electric fields. 
The only non-vanishing component of the additional electric 
field is in the first case E z = d,B y {t)x, and in the second case 
E x = -d,B y (t)z. By taking the average of the two cases, one 
obtains a more symmetric electric field, which is circular in 
the x-z plane. 

If B y (t) is constant in time, the perpendicular electric field 
vanishes and the two cases are gauge equivalent. However, 
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we do not know how to implement a constant magnetic field 
exactly in a practical numerical calculation. In the numerical 



evaluation of Eq. (63 1 one might replace the lower integration 



bound by a finite time t a . But then one effectively deals with 
a sudden switch-on perpendicular field of the following form 



By(t) 



if t < t a , 
B y if t > t a . 



(64) 



This sudden switch-on leads to large perpendicular electric 
fields, and as a result the current density in case A and B can 
be different. To instead implement a magnetic field that is 
effectively constant in time we have to switch it on slowly. In 
our numerical calculations we will take a magnetic field of the 
following form 

f if t < t a , 

B y (t) = \ B y (t - t a )/(t b - t a ) if ta < t < t b , (65) 



if t > t b . 



If we choose tb - t a large enough and % small enough before 
the important physical effects happen, the magnetic field is ef- 
fectively constant. In that situation the perpendicular electric 
fields are small and the induced current density in case A and 
B should approximately have the same magnitude. We will 
use this feature to test our methodology. 



A. A z (x) = -B y (t)x 
In the case that A z (x) = -B y (t)x we obtain 



H(t, t') = j dVx'n;; 3 



(x, x). 



(66) 



After inserting the explicit expression for S ± (x, x'), taking the 
trace, and integration over y, z, p' y and p' z we find 



H(t,t') = -q 2 Y j J] dx ' x ' 



p k'=0' 



'P 



Vp, P '(f, t')W p , P '(x,x') 



v 



(67) 



f>y=Py.Pi=Pz 



where the functions V and W are given by 

V p , p ,(r, = ImJ^W/;,^^)^,^)*, (68) 
W p , p ,(x,x') = ^gp-A^gp'sWgpAx^gp'si*')- (69) 

s=± 

The expression for H(t, t') can be simplified by performing 
the integration over x' followed by integration over p y . Using 
the relations from Appendix|A]it can be shown that 



-P- dx'x'W p , p -(x,x') 

J-oc Mi J-oo ] Py 



=Py 



2ttV2 



^k[28 Kk ,-5 k ^ k ,-8 k+uk ,}. (70) 



Inserting Eq. ( 70 1 into Eq. ( |67| i yields 

2 °° /^oo 



H(t, t') 



where 



CJ k X 



[2V Kk {t, f) - Vu-i(f, ?') - V kM1 (t, t')] p _ p , , (71) 



0) k = 



\qB z \k 



^2\qB z \k + m 2 



(72) 



To speed up the numerical computation it is convenient to take 
the wavefunctions with the same momenta together in the in- 
tegrand. For this reason we rewrite Eq. (71 1 into 

2n 2n ' 2tt f-f 2n 

[2a) k V k , k (t, t') - u) k+l V k+hk (t, t') - to k V kMl (t, t')] p =p , . (73) 



B. A x (x) = BM)z 



If A x (x) = B y (t)z we obtain 



H(t, f) = - j dV z' n 2 c l (x, x'). 



(74) 



Inserting the explicit expression for the two-point function and 
performing the trace, we find that 

Il£W) = i— 2 e i( ^^X y - y ' )+i(ft ^ )(z - z,) x 



p.p 

2 

< p ,(x,x') V^M + ^-V^itj') 

+ ^^Wl pl {x, x')V^ p ,(t, t')\ , (75) 

where we have defined the following functions 

Vi p ,{t,t') = Im^^(f)/; ! W^(0^ s (0*, (76) 

.s=± 

V^ p ,(t,t') = Im2/; J (0/^. S (0/ P + - J (f')7;- J (f')*,(77) 

.s=± 

Wp p ,(x,x') = 2j sg ps (x)gp'- s (x)g ps (x')gp,- s (x'), (78) 



p ,(x, X) - £j s gps(x)gp'-Ax)gp-s(K )gp's (x')- (79) 



To simplify H(t, t') we can make use of the following rela- 
tion 



It If 

_. r <jp z r 8 

J 2n dh 



e izh f(p z +h/2,p z -h/2) 



(80) 



/!=() 
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Furthermore, using the relations from Appendix |A| it follows 
that 



X°° dp f°° I 



: Py=Py 

\qB z \ 
In 



(<W - <W) . (81) 



r ^ r dx'w^. p ,(x,x')| , =o. 

Using the last three equations above we find 



(82) 



*=0 



+ ^ y ^i (f ' f,) ]' (83) 



where p z — p z + h/2 and p' z -p z — h/2. 



C. Numerical procedure 

We now will discuss the details of the numerical evalua- 
tion of the current density along the perpendicular magnetic 
field. Firstly, we obtain the wavefunctions fl s {t) numeri- 



for different f . 
is computed us- 



cally through solving Eqs. ( 17 i and ( 18 1 using a Runge-Kutta 
method implemented in Matlab. We make sure that we ob- 
tain fp S {t) at equally spaced time steps. The next step is to 
construct the integrand of Eqs. (73 1 and ( [83} 
The derivative in the integrand of Eq. (|83|l 
ing finite differences. We then perform the f' integration in 
Eq. (63 i using the trapezoidal rule with a lower integration 
cutoff. Thereafter we perform the p z integral, in an interval 
symmetric around p z = 0, using the trapezoidal rule. The up- 
per and lower cutoff are taken so large that varying them does 
not change the results. The last step is to perform the sum over 
k. The results are dominated by the small k values, therefore 
we sum over k until we reach convergence. Typically one only 
has to sum over a few values of k to obtain an accurate answer. 



V. RESULTS 

We will now study the current density generated by the chi- 
ral magnetic effect in parallel electric and magnetic fields with 
a perpendicular magnetic field as in Fig. 1. We will compute 
this current density numerically using the linear response re- 
lation Eq. ( |63) . Since the calculation is based on linear re- 
sponse, our results will be valid for perpendicular magnetic 
fields small compared to the parallel electric and magnetic 
fields. 

The full electromagnetic current density has two compo- 
nents. Firstly, it has a component along the perpendicular 
magnetic field due to the chiral magnetic effect. Secondly, it 
has a component in the longitudinal direction along the elec- 
tric and magnetic field due to pair production. This compo- 
nent does not vanish for weak perpendicular magnetic fields, 



and can to first order be computed using Eq. (48 1. As ex- 
plained in the introduction the longitudinal component van- 
ishes in the QCD setup due to a cancellation of the contribu- 
tion of red and green quarks. 

The electromagnetic current density will generate fields 
themselves which can modify the dynamics. This backreac- 
tion can be safely neglected as long as the fields induced by the 
currents stay small compared to the background fields. Such 
regime can always be reached by considering times short after 
the switch-on of the background fields. In this article we will 
not consider this backreaction and leave its study for future 
work. 

We will present results for the chiral magnetic effect in a 
sudden switch-on electric field and a pulsed electric field be- 
low. We consider a perpendicular magnetic field that is ef- 
fectively constant in time as in Eq. (65 1. The formalism we 



have developed in this article allows one to analyze the chiral 
magnetic effect in other settings as well. 

We have performed all numerical calculations in an effec- 
tively constant magnetic field using both gauge field choice 
A and B. These choices are approximately gauge equivalent, 
and the accuracy of the approximation can be improved by 
switching on the effectively constant magnetic field slower. 
The calculations performed using gauge field choice A and B 
are independent and hence can be used to test our methodol- 
ogy. In the numerical calculations in an effectively constant 
magnetic field we have found excellent agreement between 
the results obtained with choice A and B. 



A. Sudden switch-on electric field 

Here we will consider the chiral magnetic effect in a 
sudden-switch on electric field of the following form E(t) = 
E z 6(t). The corresponding gauge field reads A z (t) - -E z 6(f)t. 
The perpendicular magnetic field is taken to be effectively 
constant. In Ref. [ 1 1 the current density in the y-direction 
was computed exactly for t » 0. For small B y and t » the 
rate of current density generation equals iflOl 



. q 2 B y \qE z \B 2 z IB Z 



\qE z 



(84) 



In order to cancel the rapid oscillations in the current den- 
sity arising from the sudden-switch on of the electric field 
we will investigate the running average of the current density, 
here defined as 



r t+cj Vk/£J 
(Jy(t))= dt'j y (f). 

Jt-c 



t-cl \\qE z \ 



(85) 



We display the results of the numerical computation of the 
running average of the current density using c = 1 in Fig.[2]for 
different values of m = mj ^\qE z \ and B z /E z = 1. The linear 
response calculation shows that after the switch-on the current 
quickly grows linear with time. A fermion mass suppresses 



the production of chirality as can be seen from Eq. (57 1. This 
explains why the current density is smaller for particles with 
a larger mass. 
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FIG. 2: Running average of the current density generated by the chi- 
ral magnetic effect in an electric field suddenly switched on at t = 0. 
HereS z /£ z = 1. 
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FIG. 3: Rate of current density generation due to the chiral magnetic 
effect at late times in a sudden-switch on electric field as a function of 
B z /E z for different masses. The lines denote the small B y limit of the 
exact result; points indicate the results of the numerical calculations 
using linear response. 



In Fig. [3] we compare our numerical results for the rate of 
current density generation to the small B y limit of the exact 



result, given in Eq. (84 1. It can be seen that we find excellent 
agreement between the results obtained using linear response 
and the small B y limit of the exact result. Thus our linear 
response approach has passed a critical test. It implies that 
the study of the dynamics of the chiral magnetic effect using 
linear response can be performed successfully. Alternatively, 
our results can be seen as an independent verification of the 
results obtained in Ref. IflOl . 

It can be seen that the rate of current density generation in- 
creases if B z is enlarged. This is natural, since the amount of 
chirality production is increased. But at the same time enlarg- 
ing B z decreases the degree of polarization of the fermions in 
the y-direction. The combination of these two effects results 
in the saturation of the rate of current generation for large B z . 
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FIG. 4: Current density due to the chiral magnetic effect in a pulsed 
electric field, as a function of time. Here t = 1/ ^j\qE z \ and m = 0. 
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FIG. 5: Same as in Fig.|4]but now for m = 0.5. 

B. Pulsed electric field 

We will now study the chiral magnetic effect in a pulsed 
electric field that has the form E z (t) = E z / cosh 2 (f/r). The 
corresponding gauge field reads A z (t) = -E z [l + tanh(f/r)]T. 
The perpendicular magnetic field is taken to be effectively 
constant. 

In Figs. [4] and [5] we display the current density generated 
by the chiral magnetic effect as a function of time for t = 
1/ ^\qE z \ and respectively m = mj ^J\qE z \ = and m = 0.5. 
It can be seen that the current density rises quickly around 
t — 0. This is because only then an electric field is present 
so that chirality will be produced. For large t the current does 
not longer grow because there is no production of chirality 
anymore. For large values of B z /E z we find that the current 
density becomes approximately constant in time if m = 0. If 
the fermions are massive the current density exhibits a slowly 
damped sinusoidal oscillation for large values of B z /E z . The 
mass also suppresses the magnitude of the current. 

We observe that for smaller values of B z /E z the current 
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FIG. 6: Average current density at late times in a pulsed electric field 
as a function of r for different in. 



fm 




FIG. 7: The function f(m) which describes the mass dependence of 
the average current density at late times in a pulsed electric field. 
Points are numerical results, solid line is a fit. 



VI. CONCLUSIONS 

In this article we have investigated the real-time dynam- 
ics of the chiral magnetic effect using linear response theory. 
We have considered a field configuration in which a homoge- 
neous (chromo)electric field with arbitrary time-dependence 
lies parallel to a homogeneous and static (chromo)magnetic 
field. These parallel fields are the source of the chirality. To 
this field configuration we have added a perpendicular homo- 
geneous and static magnetic field. We have computed the in- 
duced current density along this perpendicular magnetic field 
explicitly for a sudden switch-on and a pulsed electric field. 

In the sudden switch-on electric field we have obtained ex- 
cellent agreement with an earlier independent analytic com- 
putation of the current density. In the pulsed electric field we 
could summarize the induced current density that we have ob- 
tained numerically with a simple analytic formula. 

The main purpose of this article was to demonstrate the dy- 
namics of the chiral magnetic effect using linear response the- 
ory. We hope that our results will be extended to other inter- 
esting field configurations in the future. For example, in heavy 
ion collisions it would be important to answer the question 
whether there is enough time for the chiral magnetic effect to 
occur in the quickly decaying magnetic field. This question 
could be addressed using our methodology. The chiral mag- 
netic effect could also be investigated with lasers that create 
strong electromagnetic fields. For that purpose it would be 
important to extend the results to field configurations that are 
as close to the experimental situation as possible. 

As a side result of our work we have obtained a derivation 
of the induced current density in static homogeneous parallel 
electric and magnetic fields. We have also obtained an ana- 
lytic formula for the chirality production for massive fermions 
in static homogeneous parallel electric and magnetic fields. 
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Appendix A: Relations involving the function g p ,(x) 

To evaluate these integrals we will make use of the follow- 
ing four identities which directly follow from the properties 
of the Hermite polynomials, 



f 



dx'x'gp^xOgp'^x') 
1 



P,=Py 



Py_ 
qB z 



^q~B z 



(yfkdksjs + y/k+s6 k+sJc >), (Al) 
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i r°° d Py 

L^ Pygp - s(x)8Mx) 



\qB z 



Py=Py 



In 



1 Hk-s6 k - 2s ,k' + ^IkhA 
^2\qB z \ 



(A2) 



f 



(Ipy 

-^gp-sWgp'sW 



Py=Py 



\qB z 
2n 



S k - Sik >. (A3) 



dxg ps (x)g p ^ s (x) 

CO 



= 8, 



(A4) 



Py=Py 



Appendix B: Wave functions in a sudden switch-on electric field 

We will review the explicit solutions for the wave functions 
fp S {t) for a sudden switch-on electric field of the form E(t) = 
E z 6(t) [22J. Then we will evaluate an integral that is necessary 
for computing the induced current and chirality production. 

For the sudden switch-on electric field we have A z (f) = 
-E z t6{t). For t < we can use the wave functions in van- 
ishing electromagnetic field, given in Eq. ( |23) . If t > the 
electric field is not longer vanishing. It follows from combin- 



ing Eqs. ( 17 1 and ( 18 i that /*(f) then satisfies 



-df-q 2 E 2 z [t + 



qE 2 



+ \sqE z 



^(0=^(0. (Bl) 



Eq. ( |B1| | is an eigenvalue equation for a particle in an upside- 
down harmonic potential. There are no bound states, so k 
is not quantized. The solution is a linear combination of 
parabolic cylinder functions D y (z), 



f*{f) = a ± s D v> {&+ftD Vi (-t), 



(B2) 



where v s = -(ssgn(qE z ) + 1 + ik 2 /\qE z \)/2 and £ = 
V2pie^ 4 (r + Pz /qE z ). 

In the conventional normalization the parabolic cylinder 
functions are given explicitly by the following integrals 



D V (Z) 



1 -e-i" I d/r-V^-r 



r(-v) 



\* r 

Jo 



(B3) 



for Re(v) < 0, and 



D v (z) = y -e? 2 ' J df f cos (zf - ±ttv) e~5' 2 (B4) 

for Re(v) > -1. Using these relations one can show that 
D v (0) = 2 V/2 V£/r[(l - v)/2]. Asymptotically the parabolic 
cylinder functions behave as follows 



1 2 V27T , 1 

lim D v (z) = z v e~* z - Cv ^^-"-i e 4 2 

|z|-t°° I (-v) 



(B5) 



with c v = for |arg(z)| < ^f,c v - exp(i7rv) for | < |arg(z)| < 
and for c v = exp(-i7rv) for — ^ < |arg(z)| < -f . These 



asymptotic functions describe the real parabolic cylinder func- 
tions very good if |z| > 2|v|. 



The constants and B* of Eq. (B2i can be found by re- 
quiring continuity of fp S (f) and its derivative at t = 0. This 
gives the following two equations 



ft = 



1 
1 

W 



Po + sp z 



Po 



[_ e W4 ^ 2 \qE z \D' Vs (-{o) 



±ip D Vs (-to)], 
+ipoA< s (&)], 



(B6) 



(B7) 



where £o = ^l2& m ^p z l\qE z \ 1 ^ 2 sgn(qE z ) and W denotes the 
Wronskian of the two independent solutions presented in 
Eq. (|B2[ ). Applying Abel's differential equation identity to 
Eq. ( |B1[ | shows that the Wronskian is independent of £ . Hence 
without loss of generality we can evaluate W at £ = 0, 
which yields W = 2e ur ^ 4 -0rjgEJ/r(— v s ). Using the relation 



D'y(z) 
to 



2zDy(z)-D v+l (z) we can simplify Eqs. (|B6|) and (|B7|) 

r(-v s ) ipo + sp z 



a, = 



V2^ 



Po 



-[A-,+i(-ft) 



Jbr/4 



y/2\qE z 



: (p z sgn( 9 £ z ) ± po) D v ,(-nft)], (B8) 



ft = 



r(-vs) /po + spz 

y/2n~ 



Po 



[A>, + l«b) 



>/4 



y/2\qE 2 



( Pz sgn(qE z )±p )Dy s (^)]. (B9) 



Using the asymptotic expansion of the parabolic cylin- 



der functions, Eq. (B4i, it is possible to obtain a very good 
approximation for the function fjj s (t) away from the points 
p z = -qE z t and p z = 0. Let us now for a moment choose 
qE z > and take y/\qE z \t » 1 . In order to obtain the induced 
current we need to evaluate \fp + (t)\ 2 , for which we find after 
taking the dominating terms in the asymptotic expansion and 
approximating p by \p z \ the following 



i/; + (or 



for p z < -qE z t - A, 



2e » + g(0 for - qE z t + A < p z < -A, 
for p z > A. 

(B10) 

where A = 2\qE z \ l ' 2 \l + ^|, f = + p 2 /<?£ z ), and 



2 1 

S(0 = -P7 



Q 4\qE z \ _ g 4|?E Z I 



xRe 



(BID 



As follows from Eq. ( |48| ) we have to evaluate the following 
integral to obtain the induced current 

/(*, = lim ^ e -«PMf [|/; +(f) |2 _ l] . (B12) 



13 



In general this integral can only be evaluated numerically. 
But using Eq. ( |B 10| > we can obtain I(k, t) exactly for large t. 
Firstly one realizes that taking the limit e — > implies that we 
have to integrate \fp + (t)\ 2 - 1 over p z in an interval symmetric 
around p z = -qE z t. The contribution to the integral in the re- 
gion where the approximation Eq. ( |B10[ ) breaks down can be 
bounded from below and above by a time-independent con- 
stant. The contributions for p z < -2qE z t and p z > will can- 
cel against each other. The non-vanishing contribution comes 
from the intermediate region in which -2qE z t < p z < 0. The 
integral over the rapidly oscillating function g(£) is subdomi- 
nant in the large t limit. The dominant contribution in the large 



t limit is qE z t coming from the interval -2qE z t < p z < -qE z t 
and qE z t[2exp(-m< 2 /\qE z \) - 1] from the interval -qE z t < 
p z < 0. Ading both contributions yields 



lim -I(K,t) 



qE z » 2 

71 



(B13) 



By repeating the analysis for qE z < one can verify that the 



sign of Eq. (B 13 1 is correct, so Eq. (B13 1 holds for all values 
of qE z . 
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